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1.0 


INTRODUCTION  AND  SUMMARY 


The  behavior  of  materials  under  particle  impact  is  of  critical  importance 
in  the  performance  of  IR  windows,  nosetips  and  heatshields  and  turbine  blade  applica- 
tions. The  applications  have  been  roughly  classed  as  high  and  low  velocity  encounters 
with  a minimum  of  interaction  between  the  techniques  used  to  investigate  performance 
in  these  different  velocity  regimes.  While  such  a separation  has  thus  far  been  prac- 
tical, the  extension  of  missile  performance  to  higher  velocities  (e.  g. , Mach  2.  5-3.  0) 
and  anticipated  performance  in  the  Mach  5.  0 and  6.  0 range,  make  it  necessary  to 
broaden  the  base  of  understanding  of  material  behavior  in  order  to  predict  material 
behavior  under  conditions  not  previously  investigated.  The  objectives  of  the  program 
are: 

1.  Provide  the  basic  understanding  of  the  particle  impact 
phenomena  which  produce  damage  in  IR  windows  in  order 
to  establish  material  development  guidelines,  predict 
performance  under  anticipated  operational  conditions 
and  define  failure  conditions. 

2.  To  establish  the  technical  base  on  which  extrapolation 
into  high  velocity  regimes  can  be  carried  out. 

3.  To  extend  the  understanding  of  material  behavior  into  the 
transitional  regime  between  low  velocity  rain  erosion 
processes  important  to  IR  and  optical  windows  and  the 
hypervelocity  reentry  conditions  which  are  critical  to 
reentry  vehicle  design. 

The  erosion/degradation  problem  can  be  and  has  been  addressed  from 
the  viewpoints  of  empirical  materials  development,  microstructural  studies,  dynamic 
code  analysis  of  the  impact  event,  and  static  (Hertzian  indentation)  analysis  and 
experimental  simulation.  Currently  lacking  is  a framework  for  logically  connecting 
these  separate  approaches,  particularly  making  the  low  to  high  velocity  transition. 
Establishment  of  such  a framework,  its  application,  and  correlation  with  experiment 
constitutes  the  basis  of  the  present  effort. 

Unique  key  elements  of  the  program  are  1)  providing  the  computational 
framework  for  translating  dynamic  code  analysis  of  the  impact  event  into  a form 


suitable  tor  mierostructurol  and  micromechanical  formulation  of  the  damage  event, 

2)  the  accompanying  microstructural  and  micromechanical  formulation  of  single  and 
multiple  impact  phenomena,  and  3)  correlation  of  the  damage  predictions  with 
existing  data. 

The  outputs  from  the  program  are  1)  improved  integration  of  the  currently 
utilized  approaches,  2)  delineation  of  the  velocity  conditions  under  which  static  linear 
formulations  become  inadequate  as  well  as  when  nonlinear  approaches  are  necessary, 

3)  coupling  of  code  solutions  to  micromechanies  and  damage  formulations,  4)  sensi- 
tivity of  damage  to  macroscopic  and  microstructural  material  parameters,  and  5)  a 
logical  formal  transition  between  single  and  multiple  impact  phenomena. 

The  benefits  from  the  program  are  1)  a baseline  technique  with  which  to 
examine  many  velocity  ranges  and  applications,  2)  a baseline  technique  capable  of 
adaptation  to  impact  code  improvements  or  other  impact  codes,  3)  an  integrated 
approach  to  formulating  material  development  guidelines,  and  4)  a formulation  suf- 
ficiently general  to  allow  low  to  high  velocity  encounter  transitions. 

Specific  tasks  pursued  under  this  contract  are  summarized  below: 

Task  I.  Impact  of  a Water  Drop  on  a Rigid  Surface 

Impact  surface  pressure  distributions  were  determined  by  the 
WAVE-L  computer  code  under  subcontract  to  California  Research  and 
Technology,  Inc.  (CaRT).  The  scope  of  the  effort  was  centered  upon 
resolving  the  early  time  off-axis  peaking  reported  by  Rochester  and 
Brunton  and  as  evidenced  in  previous  water  drop  simulations  by  CaRT. 

Task  II.  Numerical  Calculations  of  an  1100  FPS  and  4000  FPS  Impact  of 
290  Alicron  Glass  Spheres  on  a ZnSe  Plate 

The  impact  event  was  subjected  to  static  and  Hertzian  analysis. 
The  static  analysis  was  taken  from  well-known  forms  in  the  technical 
literature  and  simulated  the  theoretical  impact  pressure  distributions. 

Dynamic  analyses  were  accomplished  by  the  linear  elastic 
options  available  in  the  WAVE-L  code  which  formed  the  baseline  ana- 
lytical technique.  Two  impact  velocities  were  considered.  The  surface 
pressure  distributions  were  taken  as  those  indicated  directly  by  impact 
expe  riments. 


The  results  of  the  static  and  dynamic  computations  were 
compared  to  assess  the  range  of  validity  of  static  and  dynamic  stress 
models  and  the  significance  of  nonlinear  mechanics  to  the  velocity 
range  of  interest. 

Task  HI.  Crack  Nucleation  and  Growth  from  Low  Velocity  Particle  Impact 
on  Brittle  Materials 

Stresses  were  examined  at  various  times  by  micromechanical 
methods.  These  included  statistical  distributions  of  fracture  initiation 
sites,  multiple  crack  initiation,  slow  and  fast  crack  growth  concepts, 
and  crack  coalescence.  Micromechanical  variables  were  varied  to 
define  critical  material  concepts. 

The  outputs  of  the  analyses  were  correlated  with  available 
data  to  provide  direction  for  improved  damage  modeling,  as  well  as 
preliminary  guidelines  for  material  improvement.  Definition  of  critical 
theoretical  and  experimental  concepts  were  accomplished. 

The  following  three  sections  discuss  the  results  of  each  of  the  three  major 
tasks  of  the  Erosion  Mechanics  and  Micromechanics  Program.  Conclusions  are 
summarized  in  Section  5.0. 
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IMPACT  OF  A WATER  DROP  ON  A RIGID  SURFACE 


2. 1 Introduction 

Previous  experimental,  theoretical,  and  numerical  studies  (Refs.  1,  2, 

3)  indicate  that  when  a water  drop  impacts  a relatively  rigid  surface,  an  early  time 
pressure  peak  occurs  on  the  surface  somewhat  off  the  impact  axis.  This  pressure 
peak  may  exceed  the  "waterhammer"  pressure  (p^C^V^)  by  factors  of  two  to  three. 
This  phenomenon  is  of  practical  importance  since  the  pressure  pulse  could  activate 
crack  propagation  in  brittle  targets  by  forcing  high  pressure  water  into  surface 
openings. 

In  this  study,  a finely- re  solved  two-dimensional  numerical  solution  was 
performed  to  specifically  investigate  this  pressure  pulse  using  the  Lagrangian  finite 
difference  WAVE-L  code  (Ref.  4).  The  water  is  initially  divided  into  a number  of 
axisymmetric  volume  elements,  or  cells,  each  of  which  has  a pressure,  internal 
energy,  density,  and  particle  velocities  associated  with  it.  The  equations  for  up- 
dating the  cell  variables  during  a time  cycle  are  obtained  from  finite  difference 
analogues  of  the  differential  equations  for  conservation  of  mass,  momentum,  and 
energy,  and  the  compressible  material  properties  of  the  water.  The  material  inter- 
face between  the  water  drop  and  the  rigid  target  was  treated  as  frictionless.  A 
definitive  study  of  the  pressures  associated  with  liquid  drop  impingement  is  clearly 
desirable  and  is  the  goal  of  the  present  study. 

2.  2 Results  of  the  Numerical  Calculation 

Figure  1 shows  the  (initial)  computational  grid  used  for  describing  the 
1 mm  diameter  water  drop  impact  event.  The  impact  velocity  (Vq)  is  1100  fps. 
Figure  2 shows  the  velocity  field  (using  the  scale  bar  shown  on  the  figure)  for  mate- 
rial in  the  vicinity  of  the  impact  site  at  . 02  pi  sec  after  the  impact.  The  contact 
radius  (r^)  is  about  .008  cm  at  this  time.  Notice  the  tendency  of  the  velocity  vectors 
in  a roughly  triangular  region  to  focus  toward  the  contact  radius.  Note  also  that 
near  the  radius  of  contact  the  particle  velocity  is  not  arrested,  as  in  the  calculation 
of  the  one-dimensional  waterhammer  pressure,  but  is  deflected.  Figure  3 shows 
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Figure  3.  Velocity  Field  at  0.  03  nsec. 
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the  velocity  field  at  . 03psec.  The  focusing  of  the  velocity  vectors  toward  the 

contact  radius  is  more  evident  in  this  figure.  At  this  time  the  contact  radius  (r^) 

is  moving  outward  at  a velocity  (r  ) only  slightly  greater  than  the  sound  speed  (Cq) 

in  undisturbed  water.  As  will  be  shown  below,  the  pressure  peak  occurs  near  the 

contact  radius  and  continues  to  increase  in  magnitude  as  long  as  rc  is  greater  than 

C . Figure  4 shows  the  velocity  field  at  . 05  psec.  The  convergence  toward  the 
o 

contact  radius  is  still  qualitatively  apparent;  however,  the  influence  of  the  free 
surfaces  near  the  contact  radius  have  reduced  the  peak  pressure.  By  . 09  psec,  as 
indicated  in  Figure  5,  the  velocity  field  is  no  longer  strongly  convergent  toward  the 
contact  radius.  Figure  5 also  shows  the  very  early  stages  of  the  formation  of  a 
cylindrically  symmetric  jet.  (The  calculation  was  terminated  at  . 09  usee,  and  thus 
the  subsequent  interaction  between  the  jet  and  the  water  drop  was  not  computed.) 

The  pressure-radius  profiles  on  the  target  surface  are  shown  for  three 
times  in  Figure  6.  Note  that  the  waterhammer  pressure  (P0VQCo)  equals  5 kbars 
in  this  case.  The  pressure  near  the  axis  of  cylindrical  symmetry  remains  roughly 
at  the  waterhammer  pressure;  however,  as  time  proceeds  to  . 033  psec  a pressure 
peak  is  observed  near  the  contact  radius.  The  pressure  peak  reaches  about  2.  7 
times  the  waterhammer  pressure  at  . 033  psec. 

Figure  7 shows  pressure  versus  time  profiles  at  several  radii  on  the 
impact  surface.  The  peak  pressure  occurs  at  about  a radius  of  . 01  cm  or  .2  of  the 
droplet  radius.  Beyond  this  radius,  the  peak  stress  begins  to  fall  off.  Figure  8 
shows  the  peak  pressure  versus  radius  at  times  prior  to  . 09  psec. 

Figure  9 summarizes  the  relationships  between  contact  radius  speed 
(rc),  water  sound  speed  (Cq),  and  the  timing  of  the  maximum  pressure  predicted 
in  the  numerical  calculation.  The  rate  of  increase  of  the  contact  radius  (i*c)  as  a 
function  of  time  (t)  can  be  obtained  by  assuming  the  spherical  droplet  (radius  R)  re- 
mains undeformed  where  it  is  not  in  contact  with  the  rigid  boundary.  The  geometry 
is  indicated  in  Figure  9. 
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Figure  4.  Velocity  Field  at  0.  05  nsec. 


Figure  5.  Velocity  Field  at  0.  09  Usec. 
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Figure  7.  Pressure  Vs.  Time  at  Several  Radii  on  the  Impact  Surface 
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Figure  8.  Peak  Pressure  on  the  Impact  Surface  Vs.  Radius 


Figure  9.  Rate  of  Increase  of  Impact  Radius  and  Peak  Pressure  Vs.  Time 
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c 

The  r versus  time  curve  for  the  calculated  case  is  plotted  on  Figure  9 as  a solid 
c 

line.  Also  shown  on  this  curve  is  the  speed  of  sound  in  undisturbed  water  (Cq)  which 
is  . 15  cm/usec  or  about  5000  ft/sec.  The  dashed  curve  shows  the  peak  pressure 
versus  time  for  this  solution.  The  peak  pressure  attains  its  maximum  at  a time  of 
about  . 033  usee  which  corresponds  to  the  time  when  the  rate  of  increase  of  the  con- 
tact radius  (rc)  nearly  equals  the  speed  of  sound  in  undisturbed  water.  At  times 
greater  than  . 033  nsec,  the  peak  pressure  decreases. 

Thus,  we  conclude  that  the  peak  pressure  on  or  near  the  contact  radius 

increases  as  long  as  the  contact  radius  is  increasing  faster  than  the  wave  velocity 

in  undisturbed  water  (Cq).  Once  the  contact  radius  is  traveling  slower  than  this 

wave  speed,  then  the  pressure  waves  which  interact  with  the  water  free  surfaces 

generate  relief  waves  which  rapidly  relieve  the  "excess"  pressures  (i.  e. , pressures 

exceeding  the  1-D  shock  Hugoniot  value)  near  the  contact  radius.  By  . 05  psec,  for 

example,  the  peak  pressure  has  been  reduced  from  14  kbar  to  6 kbar.  When  the 

contact  radius  is  moving  faster  than  but  slower  than  the  local  shock  velocity 

(C  < r < C) , shock  wave  interactions  with  the  droplet  free  surface  is  still  pos- 
o c 

sible.  However,  the  compressive  stresses  generated  by  the  impact  of  droplet 
materials  on  the  rigid  surface  near  the  contact  radius  are  sufficient  to  dominate 
any  pressure  relief  from  the  shock  wave  interaction  with  the  droplet  free  surface. 

2. 3 Discussion 

The  early  experiments  of  Rochester  and  Brunton  (Refs.  1,  2)  indicated 
that  the  peak  surface  pressure  during  water  drop  impact  probably  occurred  on  the 
impact  axis  and  were  of  the  magnitude  of  0.  7 p ^C^V.  Cautions  on  this  conclusion 
because  of  the  large  size  of  the  pressure  transducer  were  noted.  Additional 
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refinement  of  experimental  technique,  including  the  development  and  use  of  smaller 

transducers  led  to  the  conclusion  that  the  pressure  of  the  center  was  still  somewhat 

less  than  the  waterhammer  pressure,  but  that  a maximum  of  1.  8 o C V occurred  at 

o o 

approximately  . 1 drop  diameter  off  the  impact  center.  In  these  experiments,  the 
sizes  of  the  pressure  transducers  were  respectively  approximately  one-fifth  and  one- 
fifteenth  of  the  drop  diameter.  The  averaging  over  transducers  of  these  sizes  con- 
tributed to  some  uncertainty  with  regards  to  the  physical  location  of  the  peak 
pressure  pulse  and  its  magnitude.  The  computed  pressure  pulse  in  Figure  (J  is  very 
sharp  spatially  and  exhibits  a maximum  whose  location  is  in  good  agreement  with 
the  more  refined  experimental  data. 

The  response  time  of  the  transducer  should  also  be  considered  in 
attempts  to  compare  the  present  calculations  with  experiment.  The  experimental 
transducers  were  estimated  to  have  a response  time  on  the  order  of  0. 1 and  0.  2 
microsecond.  As  indicated  in  Figure  7,  the  half  width  of  the  computed  pressure 
transient  during  water  drop  impact  is  somewhat  less  than  experimental  response 
time  of  the  transducer. 

Combining  the  factors  of  physical  size  of  the  pressure  transducers  and 
the  associated  response  time  leads  to  the  conclusion  that  the  calculated  maximum 
of  2.  7 PQC0V0  is  not  in  disagreement  with  experimental  limitations.  Further  sub- 
stantiation of  this  experimental  discretization  effect  w'as  evidenced  in  earlier  cal- 
culations performed  with  a coarser  computational  mesh.  Careful  interrogation  of 
those  results  showed  a pressure  maximum  occurring  off-axis  in  a similar  fashion. 

The  computed  magnitude  of  this  pressure  pulse  was  lower,  about  1.  6 p C V,  and 

o o 

there  was  then  uncertainty  as  to  the  validity  of  such  narrow  transients. 

Examination  of  the  pressure  time  distribution  given  in  Figure  7 for  the 
center  of  the  impact  leads  to  the  impression  that  the  present  finely  zoned  calculations 
predict  greater  pressures  than  experimentally  reported.  However,  if  the  relatively 
long  experimental  rise  times  for  pressure  transducers  (»■  0.  1 psec)  are  taken  into 
account,  the  appropriate  time  average  of  the  present  computations  would  not  be 
inconsistent  with  the  experimental  observation. 
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The  present  studies  provide  guidelines  for  the  formulation  of  experi- 
mental programs  to  refine  the  understanding  of  the  impact  pressures.  If  the  present 
computations  are  representative  of  actual  impact,  then  pressure  transducers  with 
fine  spatial  and  temporal  resolution  are  required.  The  discretization  similarities 
between  coarsely  zoned  computations  and  transducers  of  finite  size  should  be 
recognized  in  experimental  planning  and  interpretation  of  data.  The  results  given 
in  Figures  6 and  7 indicate  that  the  transducers  should  be  quite  small  (only  a very 
few  percent  of  the  drop  diameter)  and  the  combined  rise  time  of  the  transducer  and 
instrumentation  should  be  less  than  0.  01  nsec.  The  effects  of  spatial  and  temporal 
resolution  could  be  simulated  by  time  and  space  averaging  of  the  present  results, 
but  such  a process  is  outside  the  scope  of  the  present  study. 

Heymann  has  cited  several  experimental  results  (Refs.  5,  6,  7)  which 
would  be  consistent  with  his,  and  the  present,  prediction  of  a ring  of  off-axis  pres- 
sure in  excess  of  the  one-dimensional  waterhammer  pressure.  In  common  with  the 
findings  of  Heymann  and  Rochester  and  Brunton,  the  zone  of  high  impact  pressure 
may  be  described  as  "localized  and  fleeting." 

The  relationship  between  the  present  pressure  computations  and  localized 
material  damage  is  addressed  later.  This  division  is  deliberate  as  the  physical 
principles  governing  material  response  such  as  crack  or  failure  initiation  and  growth 
need  not  be  invoked  to  describe  the  present  results.  What  is  of  importance  here  is 
that  fluid  flow'  characteristics  account  for  a pressure  build-up  phenomena  of  a differ- 
ent character  than  that  associated  with  low  speed  solid  impact,  and  may  provide  the 
mechanism  of  differing  materials  response  between  liquid  and  solid  impact.  Penetra- 
tion of  fluid  flow  into  minute  surface  cracks  would  certainly  be  sensitive  to  the 
impact  pressure  distribution.  This  mechanism,  for  example,  would  have  to  be 
considered  in  evaluating  materials  response  to  the  pressures  developed  during  hydro- 
impact. 
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2.4 


Summary 


The  present  computations  represent  the  most  finely  zoned  hydrodynamic 
computation  yet  reported  for  liquid  drop  impact.  These  computations  were  specif- 
ically tailored  to  resolve  the  experimental  and  theoretical  uncertainties  associated 
with  pressure  distributions  developed  during  hydrodynamic  impingement.  This 
objective  was  believed  accomplished  as  evidenced  by  the  well  behaved  flow  fields 
obtained  and  the  easy  physical  explanation  of  the  computational  results.  It  is  con- 
cluded that  significant  highly  transient  pressure  pulses  exist  off  the  impact  center 
in  general  accordance  with  Rochester  and  Brunton's  (Ref.  1,  2)  experiments.  The 
pressure  peak  occurs  when  the  velocity  of  the  radius  of  contact  approximates  the 
speed  of  sound  in  water.  The  basic  nature  of  the  phenomena  appears  general,  and 
should  be  considered  as  a possible  basis  for  differences  between  material  damage 
mechanisms  and/or  observations  for  liquid  and  solid  impact. 


3.0  NUMERICAL  CALCULATIONS  OF  AN  1100  FPS  AND  4000  FPS 
IMPACT  OF  290  MICRON  GLASS  SPHERES  ON  A ZnSe  PLATE 

3. 1 Introduction 

The  numerical  analysis  of  the  dynamic  velocity  and  stress  fields 
resulting  from  the  impact  of  glass  spheres  impacting  elastic  ZnSe  was  performed 
using  a Lagragian  finite  difference  code.  A prime  objective  was  to  compare  the 
dynamic  results  of  this  numerical  analysis  with  the  quasi-static  results  predicted 
by  a Hertzian  analysis.  The  purpose  is  to  indicate  the  nature  and  significance  of 
departures  from  the  quasi-static  model  concepts  commonly  used  in  the  literature. 

, The  elastic  pr-ojectile  and  target  are  divided  into  a number  of  volume 

elements,  or  cells,  each  of  which  has  a set  of  stresses,  velocities,  etc. , associ- 
ated with  it.  The  equations  for  updating  the  cell  variables  during  a time  cycle  are 
obtained  from  the  finite  difference  analogues  of  the  conservation  of  mass,  momentum, 
and  energy  equations  along  with  the  elastic  constitutive  relations  for  the  materials 
involved.  Table  1 indicates  the  elastic  properties  used  in  this  study. 

TABLE  1.  ELASTIC  PROPERTIES  OF  GLASS  SPHERE 
AND  ZnSe  TARGET 


Density 

Poisson's 

Young's 

PQ  (gm/cc) 

V 

E (Mbar) 

Glass  Projectile 

2.23 

.2 

.64 

ZnSe  Target 

\ 

5.27 

.3 

.69 

\ 
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3.2 


Results 


f* 

Figure  10  shows  the  initaal  computational  grid  used  for  both  the  1100 
and  1000  fps  cases.  Figures  11,  12  and  13  show  principal  stresses  (magnitude  and 
direction)  for  the  1100  fps  impact  case  at  a time  of  . 14  psec.  The  plots  roughly 
correspond  to  the  time  of  maximum  penetration.  The  Hertzian  theory  predicts  a 
time  of  peak  penetration  of  . 15  psec,  while  the  dynamic  solution  predicted 
. 13  usee. 


c*ax  ten 


Radius  (cm) 


Figure  10.  Computational  Grid  Near  the  Impact  Point 
for  the  1100  and  1000  fps  Impacts  of  a 
290  - . 029  cm  Glass  Sphere  on  a ZnZe 

Target 
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Figure  11.  Comparison  of  the  Stress  Fields  for  the  Dynamic  Solution  (V  - 1100  fps), 
: nd  for  the  Hertzian  Approximation  at  . 14  psec  (T/T  = 1.  0!j^ 


Figure  12.  Comparison  oi  the  In- Plane  Tensile  Stresses  for  the  Dynamic 
Solution  (VQ  - 1100  ips)  and  for  the  Hertzian  Approximation  at 
. 14  pi  sec  (T/T  =1.05) 


Figure  13.  Comparison  of  the  Out-of-Plane  Tensile  Stresses  for 
the  Dynamic  Solution  (1100  fps)  and  for  the  Hertzian 
Approximation  at  . 14  psec  (T/T  1.  05) 


Figure  11  is  a comparison  of  the  Hertzian  stress  field  and  the  calculated 
stress  field.  The  compressive  stress  fields  are  similar  except  near  the  projectile. 

As  discussed  below,  the  radial  and  hoop  stresses  are  much  lower  in  the  dynamic 
calculation  near  the  axis  of  symmetry.  Note  that  there  are  no  stresses  in  the 
projectile  in  the  Hertzian  case,  whereas  the  dynamic  calculation  included  the  elastic 
projectile  and  the  stress  reflections  and  interactions  between  the  projectile  and 
deforming  target  surface.  Figure  12  compares  the  tensile  stresses  in  the  r-z  plane 
as  predicted  in  the  Hertzian  quasi-static  model  with  the  numerically  calculated  in- 
plane tensile  stresses.  The  dynamic  tensile  stresses  are  higher  in  magnitude  and 
extend  closer  to  the  projectile  as  compared  to  the  Hertzian  case.  Figure  13  com- 
pares the  out-of-plane  or  hoop  tensile  stresses  in  the  Hertzian  versus  dynamic  case. 
The  hoop  tensile  stresses  are  again  larger  and  again  extends  closer  to  the  penetrator. 
Also,  note  that  hoop  tensions  exist  in  the  dynamic  case  at  the  target  front  surface 
about  two  radii  from  the  axis  of  symmetry.  In  the  Hertzian  quasi-static  solution,  no 
hoop  tensions  develop  near  the  surface.  We  therefore  conclude  that  the  dynamic 
solution  represents  a much  more  severe  tensile  environment  than  does  the  Hertzian 
solution. 

Figures  14,  15,  and  1G  show  comparisons  for  the  1100  fps  case  of  the 
components  of  the  stress  tensor  near  the  target  surface  at  three  times.  Figures  17 
through  22  show  the  same  stress  comparisons  at  the  time  of  peak  penetration  for  the 
4000  fps  case.  Also  on  these  plots  are  a comparison  with  the  predictions  of  the 
Hertzian  theory.  The  agreement  with  respect  to  the  axial  stress  (a  ) is  good.  This 
agreement  is  responsible  for  the  good  comparison  of  time  to  peak  penetration  depth 
and  actual  penetration  between  the  calculations  and  Hertzian  theory’.  Table  2 summa- 
rizes these  comparisons.  The  radial  and  hoop  components  of  stress  and  a^g), 
however,  have  a distinctly  different  character  in  the  dynamic  solution  as  compared  to 
the  Hertzian  theory.  The  differences  concern  the  behavior  directly  under  the  pro- 
jectile. The  Hertzian  theory  predicts  a peak  compressive  radial  and  hoop  stress  near 
the  axis,  while  the  dynamic  solution  shows  the  radial  and  hoop  stresses  approaching 
zero  stress  near  the  axis  of  cylindrical  symmetry.  This  difference  is  due  to  the 
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Figure  1G.  Stress  Field  Along  the  Surface  of  the  ZnSe  Plate  at  . 22  (jse 
(T/Tp  - l.G)  for  the  Dynamic  Solution  (1100  fps)  and  for  a 
Hertzian  Impact 


Figure  17.  Velocity  Field  and  Computational  Guide  for  4000  fps  Solution  at  .07 
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Figure  18.  Stress  Field  in  ZnSe  Plate  for  4000  fps  Less 
Hertzian  Stress  Field  at  .07  psec 
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Figure  20.  Velocity  Field  and  Computational  Grid  for  4000  fps  Solution  at  . 11  |asec 
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Figure  21.  Stress  Field  in  ZnSe  Plate  for  '1000  fps  Solution 
Less  Hertzian  Stress  Field  at  . 11  jjsec 
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Figure  22.  Comparison  of  the  Stress  Fields  for  4000  fps  Solution  and 
for  the  Hertzian  Approximation  at  . 11  jjsec 


TABLE  2.  COMPARISON  BETWEEN  DYNAMIC  SOLUTION 
AND  HERTZIAN  THEORY 


Dynamic  Solution 

Quasi-static 
Hertzian  Theory 

Peak  Normalized  Contact  Radius 

A /R 

max 

.55 

. 49 

Peak  Normalized  Penetration, 

a /R 

max 

.21 

.24 

Time  of  Peak  Contact  Radius, 
T 

P 

. 13  nsec 

. 15  jjsec 

deformation  of  the  target  front  surface  caused  by  the  spherical  penetrator.  This 
deformation  is  realistically  treated  in  the  dynamic  calculations,  but  ignored  in  the 
small  displacement  Hertzian  theory. 

3. 3 Discussion 

The  calculated  results  of  a Mach  1 impact  of  a 290u  glass  bead  on  a 
ZnSe  target  were  examined  in  detail.  The  impact  conditions  were  chosen  to  comple- 
ment previous  studies  performed  by  Ca.  R.T.  for  AFML.  The  major  differences 
are  the  diameter  of  the  impacting  particle,  and  the  suppression  of  the  fracture  con- 
ditions within  the  target.  In  other  words,  this  is  the  dynamic  equivalent  of  the  Hertz 
impact  simulation. 

There  are  several  levels  of  comparison  to  consider.  First,  consider 
the  indentation  parameters  of  maximum  penetration,  radius  of  contact,  and  arrest 
time  (or  time  of  maximum  penetration)  for  the  particle  (Table  2).  Nothing  blatently 
different  is  evident.  In  fact,  from  this  point  of  view,  one  is  tempted  to  conclude  that 
at  Mach  1,  the  Hertz  model  does  quite  well.  Tins  has  been  found  to  be  a premature 
conclusion. 


Tlie  differences  in  Table  2 can  be  very  easily  rationalized.  The  larger 
contact  radius  of  the  dynamic  problem  is  to  be  expected  because  of  the  large  defor- 
mations at  the  point  of  impact.  As  is  true  of  many  linear  elasticity  problems,  the 
Hertz  boundary  conditions  are  satisfied  with  the  implicit  notion  of  small  displace- 


ments. A sphere  indented  to  about  one-fourth  of  its  radius  results  in  large  defor- 


mations of  the  target  and  sphere  so  one  should  expect  the  contact  radius  to  be  larger 
than  that  calculated  for  the  Hertz  model.  Because  of  this  larger  pressure  bearing 
area,  one  should  expect  the  maximum  penetration  calculated  by  the  dynamic  WAVE-L 
method  to  be  somewhat  smaller.  Of  course,  granted  the  above  two  considerations, 
the  natural  conclusion  is  that  the  time  needed  to  bring  the  impacting  sphere  to  rest 
should  be  somew'hat  shorter  because  of  the  smaller  penetration.  All  of  the  above  are 
consistent  with  the  WAVE-L  versus  Hertzian  comparison. 


The  comparison  of  the  surface  pressures  (Figure  23)  reveals  additional 
reasonable  correlation  consistent  with  the  above.  So  far,  then,  the  Hertz  simulation 
does  not  look  bad  even  for  Mach  1 impact. 


If  the  principal  stresses  are  compared  for  times  approximately  half  way 
to  maximum  and  maximum  indentation  (Figures  24  and  25),  not  much  strikes  the  eye. 
Closer  inspection  reveals  that  excepting  the  compressive  wave  propagation  phenomena, 
the  hoop  or  circumfei'ential  direction  stresses  are  the  major  visible  differences. 


These  comparisons  do  not,  however,  tell  the  correct  story  as  to  the 
behavior  of  matei'ials  under  impact  conditions.  The  figures  are  scaled  by  the  com- 
pressive stresses  which  are  much  larger  than  any  tensions  developed.  The  material 
response  of  cracking  and  degradation  of  transmission  is  dominated  by  the  in-plane 
tensions,  and  the  tension  comparison  is  the  correct  one  for  evaluating  the  applicability 
of  the  Hertz  model,  not  the  compression  dominated  indentation  parameters. 


'Hie  comparison  is  striking  (Figures  2G  and  27).  In  these  figures,  com- 
pression is  suppressed  from  plotting  and  stress  scaling  for  explicit  comparison  of 
tensions  has  been  performed.  Not  a lot  of  differences  are  noted  for  stresses  developed 
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Figure  23.  Stress  Normal  to  Indentor  Versus  Radius  at  0.  037  usee 
for  Dynamic  Solution  and  Accroding  to  Ilertizian  Theory 
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at  the  surface.  However,  very  large  differences  are  noted  in  the  interior  of  the  target. 
Probably  the  most  significant  aspect  of  this  is  that  even  if  fracture  initiation  of  the 
free  surface  could  be  suppressed  by  surface  treatments  such  as  etching  or  even  thin 
coatings  sized  taking  into  account  the  rapid  fall-off  of  radial  tensile  stresses  with 
depth  intrinsic  to  the  static  Hertz  model,  the  subsurface  tensions  which  arc  revealed 
in  the  dynamic  problem  will  contribute  to  cracking  and  degradation  of  optical  response 
if  not  outright  fracture. 

Before  accepting  this  conclusion,  however,  one  should  capsulize  the 
reasons  for  the  dynamic  versus  static  analysis  disparity.  The  first  clue  is  its  rela- 
tion to  a wave  which  moves  slower  than  the  compressive  (i.e. , "p")  wave,  in  other 
words  the  shear  or  "s"  wave.  Shear  is,  after  all,  identical  to  equal  and  mutually 
perpendicular  tension  and  compression. 

The  simple  diagram  in  Figure  28  extracts  the  basic  mechanical  phenomena 
of  interest.  The  axisymmetric  pressure  distribution  is  replaced  by  symmetric  (axi- 
svm metric  if  you  like)  point  loads  (step  loading  with  time  to  be  more  precise).  The 
p and  s waves  propagate  as  showai  and  the  sign  of  the  shears  developed  is  shown. 

The  situation  at  a later  time  is  indicated  in  Figure  29.  An  essential  feature  is  that 
as  the  shear  waves  approach  the  axis  of  symmetry,  they  have  opposite  signs.  On 
the  axis  of  symmetry',  they  will  superimpose  and  cancel.  Off  the  axis  of  symmetry, 
zones  can  be  defined  in  which  no  superposition  has  occurred  (i.e. , only  one  "s"  wave 
has  passed  the  point)  as  well  as  regions  in  which  superposition  but  incomplete  cancel- 
lation has  occurred  (i.  e. , points  through  which  both  "s"  waves  have  passed  but  because 
they  have  traveled  unequal  distances  from  their  origins,  the  wave  front  strengths  are 
unequal  due  to  inverse  powerlaw  attenuation  effects). 

The  result  of  this  wave  front  propagation  is  that  in  general  there  will  be 
zones  in  which  shear  stresses  do  not  reflect  the  static  symmetry,  and  are  higher  than 
the  static  solution.  Since  shear  is  a combination  of  tension  and  compression,  the 
dynamic  in-plane  tensions  arc  substantially  higher  than  the  static  solution. 
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Figure  29.  Superposition  of  Shear  Fronts 
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CRACK  NUCLEATION  AND  GROWTH  FROM  LOW  VELOCITY 
PARTICLE  IMPACT  ON  BRITTLE  MATERIALS 


4.  1 Introduction 

For  low  velocity  impact  of  solid  spherical  particles,  the  quasi-static 
Hertz  model  well  describes  the  stress  state  and  deformation  associated  with  maxi- 
mum indentation.  There  are  general  correlations  between  the  Hertzian  tensile 
stresses  for  a variety  of  impact  conditions,  the  material  tensile  strength  and 
tensile  strength  statistics,  and  the  observation  of  cracking  and/or  mass  loss. 

There  are  also  a variety  of  correlations  between  fatigue  behavior  and  the  multiple 
impact  mass  loss  and  cracking  behavior. 

The  purpose  of  this  task  is  to  present  a basis  for  the  prediction  of  the 
combined  effects  of  crack  nucleation  and  growth,  and  to  apply  this  formalism  to 
liquid  and  solid  impact  data.  The  basic  working  hypothesis  is  that  the  probability 
of  observing  a crack  as  a result  of  impact  is  equal  to  the  joint  probability  that  a 
crack  was  initiated  (or  a crack  source  was  activated)  and  the  probability  of  growth 
to  a size  or  geometry  which  is  experimentally  detectable  (e.g. , direct  microscopy 
or  fracture  strength  reduction). 

4. 2 Background 

The  two  basic  concepts  to  be  used  are  probabilities  of  activation  of 
crack  sources  and  probabilities  of  crack  growth.  These  will  be  illustrated  in  the 
following  discussions. 

Activation  of  crack  sources  was,  perhaps,  first  discussed  by  Griffith  in 
terms  of  the  unstable  propagation  of  pre-existing  cracks.  Basically,  this  work  con- 
sisted of  careful  delineation  of  the  crack  tip  stress  state  and/or  gross  energy 
balances  for  the  system  as  a function  of  crack  size  and  applied  stress.  The  condi- 
tions for  crack  growth  were  definable  as  either  the  exceedance  of  a local  crack  tip 
stress  state,  or  of  a crack  tip  strain  energy  release  rate.  The  instability  resulted 
from  the  fact  that  once  this  condition  was  met  or  exceeded,  further  crack  growth 
would  guarantee  that  the  condition  would  still  be  exceeded.  Time  was  not  an  explicit 
field  variable. 


Weibull  applied  the  concept  of  "weakest-link"  to  the  fracture  of  solids 
to  describe  the  mean  value  strength,  the  strength  statistics,  and  strength  increases 
for  small  stressed  regions.  This  concept  is  more  of  a mathematical  statistical 
model  than  a physical  one.  The  weak  links  are  described  by  convenient  mathematical 
distribution  functions  and  parameters  which  allow  the  fitting  of  statistics  of  fracture 
(as  well  as  a host  of  other  data  not  involving  the  weakest  link  postulate). 

Fisher  and  Holman  introduced  the  concept  of  flaw  size  and  orientation 
distributions,  which  when  coupled  to  the  Griffith  relationship,  provided  a statistical 
failure  criterion  explicitly  coupled  to  flaw  characteristics.  This  approach  also  pre- 
dicted a strength  versus  size  relationship  through  the  weak-link  postulate  that  the 
failure  probability  was  equal  to  the  probability  of  propagation  of  one  or  more  cracks, 
rn  common  with  the  Griffith  and  Weibull  approaches,  explicit  considerations  of  crack 
growth  kinetics  were  taken  into  account. 

The  early  studies  of  crack  growth  related  to  impact  phenomena  were 
associated  with  attempts  to  define  the  "equations  of  motion"  for  a crack  in  a con- 
tinuum. In  addition,  studies  have  been  conducted  from  an  atomistic  or  statistical 
mechanics  point  of  view.  Some  of  these  methods  also  introduce  explicit  considerations 
of  microstructure. 

It  is  probably  fair  to  comment  that  the  problem  of  the  kinetics  of  crack 
propagation  is  very  complex,  and  that  much  additional  work  is  required  in  this  area. 
Fatigue  behavior  is  probably  in  a similar  status.  Short  term  objectives  for  progress 
in  these  areas  arc  the  development  and  application  of  empirical  correlations  of 
laboratory  studies. 

Some  success  has  been  attained  in  correlating  crack  propagation  studies 
with  the  stress  intensity  concept;  i.  e. , the  character  of  the  elastic  stress  state 
near  the  crack  tip.  For  slow  crack  propagation,  the  following  empirical  relation- 
ships have  been  found  useful. 

clii 

= f(a/a)  = f(k) 
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where 

k = c/a  = index  of  crack  tip  stress  intensity 
a = applied  stress  field 
a = crack  size 

and  common  forms  are 


or 


da 

dt 


A (a/a) 


da 

dt 


B exp 


c/a 

KT 


The  above  correlations,  of  course,  are  limited  to  experimental  geometries  where 
the  stress  intensity  factors  are  well  known,  and  generally  to  crack  velocities  con- 
siderably less  than  the  elastic  wave  speeds  in  the  material  under  consideration. 

Tn  many  nonmetallic  materials,  the  range  of  stress  intensities  (i.e. , combinations 
of  crack  size  and  stress)  over  which  moderate  velocity  crack  growth  occurs  is  very 
narrow. 

For  such  materials,  it  is  convenient  to  define  an  approximate  critical 
stress  intensity  above  which  unstable  crack  growth  occurs  at  a high  velocity  and 
below  which  essentially  no  growth  will  occur.  In  an  analogous  manner,  some 
materials  exhibit  a small  range  of  breaking  strengths,  and  for  such  materials  it  is 
convenient  to  define  a single  material  strength. 

In  the  present  work,  the  intent  is  specifically  not  to  make  either  of  these 
two  assumptions.  The  consequence  is  to  introduce  a certain  amount  of  additional 
mathematical  labor,  but  allow  the  description  of  additional  aspects  of  the  impact 
process.  The  benefits  and  shortcomings  of  this  approach  is  discussed. 

4.  3 Method  of  Micromechanics  Analysis 

4.  3.  1 Stress  Analysis  of  Impact 

As  discussed  previously,  the  Hertz  model  provides  a good  description  of 
the  peak  stress  distribution  associated  with  solid  spherical  impact.  This  approximation 
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is  a good  model  if  the  ratio  of  the  contact  radius  to  contact  time  is  small  relative 
to  the  shear  and  longitudinal  wave  velocities.  Low  velocity  impact  will,  therefore, 
be  assumed  to  be  well  described  by  the  Hertz  model  for  the  purposes  of  this  study. 

4.3.2  Crack  Propagation 

The  crack  growth  will  be  approximated  by  the  empirical  crack  velocity 
versus  crack  tip  stress  intensity  relationships.  For  con'-enience,  the  polynomial 
form  is  used.  This  specific  relationship  has  the  disadvantage  that  imder  constant 
stress  conditions  for  a finite  time,  infinite  crack  growth  rates  and  sizes  can  be 
predicted.  To  circumvent  this  mathematical  characteristic,  the  crack  propagation 
velocity  is  taken  as  the  shear  wave  speed  if  the  empirical  relationship  predicts  crack 
velocities  in  excess  of  the  shear  wave  speed. 

This  modification  is,  admittedly,  somewhat  arbitrary,  but  is  consistent 
with  the  concept  of  cracks  rapidly  accelerating  to  a speed  approximating  the  shear 
wave  speed.  An  alternative  approach  of  fitting  the  empirical  data  to  a mathematical 
relationship  which  contained  an  upper  limiting  crack  velocity  would  produce  quali- 
tatively similar  results.  With  the  current  state  of  knowledge  on  the  theory  of  crack 
growth,  arbitrariness  must  be  introduced  in  any  attempt  to  extrapolate  observed 
low  velocity'  to  high  velocities.  However,  because  of  the  rapid  transition  to  sonic 
velocities  of  propagation,  such  details  may  assume  lesser  importance  than  implied 
by  first  consideration. 

4.  3.  3 Crack  Initiation 

Three  approaches  may  be  taken  for  this  aspect  of  the  problem.  The  first 
is  to  assume  the  material  is  free  of  microstructural  features  such  as  microcracks, 
and  that  microcracks  are  nucleated  by  mechanisms  such  as  crystallographic  slip 
in  grains,  stress  concentrations  around  pores,  inclusions,  or  agglomeration  of  such 
microstructural  features.  The  second  is  to  assume  that  the  effect  of  such  nucleation 
processes  is  an  initial  crack  size.  The  third  is  to  assume  that  the  cracks  are  pre- 
existing, albeit  very  small.  For  many  purposes,  the  second  and  third  approaches 
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are  equivalent.  The  nucleation  approach  has  appeal  for  formulating  materials  of 
complex  microstructure,  but  the  prediction  of  size  of  the  nucleated  crack  is  not 
straightforward. 

The  approach  taken  is  that  the  cracks  are  treated  as  pre-existing,  and 
that  the  size  distribution  is  related  to  the  failure  strength  distribution  observed  in 
slow  loading  environments;  e.g. , tensile  testing  and/or  static  indentation.  The 
specific  details  are  considered  to  be  obtained  from  the  integration  of  the  crack  growth 
relationships.  The  condition  of  fracture  for  a uniformly  stressed  material  is  that 
the  crack  should  grow  to  a size  comparable  with  the  dimensions  of  the  specimen. 

For  the  power  law'  of  crack  velocity,  it  is  convenient  to  define  fracture  as  equivalent 
to  predicted  crack  growth  to  infinite  size.  Under  constant  loading  rate  conditions, 
the  fracture  condition  may  be  written  as 


For  large  values  of  the  exponent  in  the  crack  velocity  equations  (i.  e. , 
sharp  transition  between  slow  and  rapid  growth),  the  strength  is  approximately  pro- 
portional to  the  inverse  square  root  of  initial  crack  size,  and  is  not  highly  dependent 
upon  the  loading  rate.  The  Griffith  criterion  is  well  satisfied  under  these  conditions. 
As  the  exponent  n decreases,  the  transition  from  slow  to  rapid  growth  occurs  over 
a wider  range  of  test  conditions,  the  dependence  of  strength  on  load  rate  becomes 
more  significant,  and  the  Griffith  criterion  becomes  a poorer  approximation.  Quali- 
tatively, similar  considerations  will  result  from  other  crack  velocity  relationships. 
The  essential  point  of  the  discussion  is  that  if  the  dynamics  of  crack  gi*owth  are 
known,  then  the  strength  and  strength  distributions  of  the  material  can  be  used  to 
infer  characteristics  of  the  flaw  size  distribution. 
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t.  4 Micromechanics  Analyses 

4.4.1  Solid  Versus  Liquid  Particle  Impact 

The  stress  analyses  have  confirmed  the  Rochester  and  Brunton  experi- 
ments (Refs.  1,  2)  which  indicated  that  off-axis  pressure  peaks  occur  as  a result  of 
liquid  (not  solid)  particle  impact.  The  phenomenon  is  apparently  real  and  there  is 
a simple  explantion  (Section  2.  0).  Therefore,  the  mechanics  of  liquid  and  solid 
impact  are  quantitatively  different  which  is  an  extremely  important  conclusion  for 
guiding  material  test  development.  This  can  be  built  upon  through  micromechanical 
considerations  to  relate  the  observed  differences  in  material  response. 

The  mechanics  were  aptly  summarized  in  Figure  9,  wherein  the  rate  of 
increase  of  contact  radius,  speed  of  sound  in  water,  and  peak  surface  contact  pres- 
sure are  given  as  a function  of  time.  The  peak  off-axis  pressure  occurs  when  the 
speed  of  sound  and  the  velocity  of  surface  contact  are  approximately  identical.  Under 
these  conditions,  the  physical  events  occurring  at  the  outer  contact  point  are  sensi- 
tive to  the  pressures  developed  earlier  (more  central  to  the  contact  area)  and  a pres- 
sure superposition  occurs  which  gives  rise  to  the  exceedance  of  the  water  hammer 
pressure. 

The  relationship  to  materials  response  is  as  follows.  Building  on  the 
concept  that  in  liquid  impact,  the  surface  pressure  drives  the  propagation  of  surface 
crack  nuclei.  Figure  30,  the  concept  that  crack  growth  obeys  laws  similar  to  those 
found  for  crack  growth  studies.  Figure  31,  the  idea  that  the  driving  force  for  solid 
impact  crack  growth  is  the  surface  Hertizian  stress,  Figure  32,  then  the  solid  and 
liquid  impact  damage  initiation  (onset  of  mass  loss)  by  Adler,  Figure  33,  can  be 
reconciled.  This  reconciliation  is  given  in  Figure  34. 

Figure  34  shows  the  improved  theoretical  correlation  with  the  solid  and 
liquid  impact  experiments  which  resulted  from  the  analyses. 

While  this  is  not  a correlation  for  all  of  the  observations  on  similari- 
ties and  differences  between  solid  and  liquid  impact,  more  detailed  coupling  of  the 
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mechanics  and  micromechanics  will  provide  the  required  understanding.  Much  of 
this  methodology  can  result  from  further  interrogation  of  the  present  computations. 

4.4.2  Surface  Finish 

L.  A.  Glenn  (Ref.  8)  experimentally  examined  the  impact  of  steel 
cylinders  on  SIA1  and  20  percent  PbO  glass  blocks  with  and  without  surface  etching 
to  maiify  brittle  fracture  strength  200 p material  removal).  Impacting  conditions 
were  14  mm  diameter  x 14  mm  long  cylinders  at  100  and  300  meters/sccond.  The 
fracture  patterns  for  etched  versus  unetched  were  quite  different  in  shape  and  size. 
Asymmetric  fracture  was  common,  but  is  probably  related  to  the  experimental  dif- 
ficulty of  getting  good  impact  alignment  with  cylindrical  projectiles.  There  is  a 
strong  possibility  that  it  has  some  association  with  the  flaw  statistics.  The  etched 
fractures  are  related  to  fracture  initiation  in  depth  and  the  unetched  fractures  are 
related  to  surface  nucleated  fracture. 

There  are  many  implications  stemming  from  the  present  comparison. 

It  is  well  known  by  the  producers  of  IR  domes  that  surface  finish  is  one  of  the  more 
important  aspects  of  production  and  by  some  small  companies  is  quite  proprietary. 
The  good  surface  finish  can  cost  more  than  basic  dome  fabrication.  Another  appli- 
cation is  that  as  impact  velocity  is  increased  beyond  Mach  1,  the  dynamical  (wave 
propagation)  phenomena  will  become  more  and  more  evident,  and  the  deficiencies 
of  the  Hertz  model  or  any  static  model  as  a basis  for  data  correlation  will  become 
more  prominent. 
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CONCLUSIONS 


The  following  conclusions  have  been  reached  as  a result  of  the  Erosion 
Mechanics  and  Micromechanics  Program: 

• Dynamic  stress  analyses  performed  under  this  contract 
have  demonstrated  correlation  with  Rochester  and  Brunton's 
experiments  which  indicated  that  off-axis  pressure  peaks 
occur  during  liquid  particle  impact.  Therefore,  the  mech- 
anics of  impact  are  significantly  different  for  liquid  versus 
solid  particles  and  can  be  expected  to  influence  material 
response. 

• At  low  velocity,  solid  particle  impact  varies  from  liquid 
particle  impact  in  the  following  ways: 

• Solid  particles  produce  higher  and  more  localized 
peak  surface  pressures  on-axis. 

• Solid  particles  produce  higher  magnitudes  of 
pressure  acting  over  shorter  time  periods. 

• Solid  particles  produce  higher  tensile  stresses 
acting  over  shorter  times. 

• Material  performance  ranking  must  account  for  the  coupling 
of  mechanics  with  material  response  at  both  macroscopic 
and  microscopic  levels  and  requires  a coordinated  effort 
which  can  treat  the  various  aspects  in  a logical  self-consistent 
fashion. 

• Surface  cracks  effect  the  material  response  during  water  drop 
impact.  Therefore,  surface  finish  is  an  important  consideration. 

• Techniques  are  available  for  merging  macromechanics  (impact 
stress  and  strain  computations)  to  micromechanics  (materials 
response,  crack  initiation,  crack  propagation)  to  provide 
improved  material  development  guidelines. 

• The  crack  propagation  considerations  provide  the  fundamental 
relationship  between  fatigue  behavior  and  long  term  erosion 
behavior. 

• Surface  finish  can  provide  enhanced  fatigue  resistance. 
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• The  best  approach  for  crack  initiation  modeling  is  to 
treat  the  cracks  as  pre-existing  and  to  infer  their  size 
and  orientation  distributions  (statistical)  from  the 
failure  strength  distributions  measured  in  slow  loading 
environments. 

• In  water  impact,  the  peak  off-axis  pressure  occurs 
when  the  speed  of  sound  in  the  water  equals  the  velocity 
of  surface  contact  radius. 
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